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Abstract 

We study the behavior of Lagrangian perturbative solutions. For a spherical void model, the 
higher order the Lagrangian perturbation we consider, the worse the approximation becomes in late- 
time evolution. In particular, if we stop to improve until an even order is reached, the perturbative 
solution describes the contraction of the void. To solve this problem, we consider improving 
the perturbative solution using Shanks transformation, which accelerates the convergence of the 
sequence. After the transformation, we find that the accuracy of higher-order perturbation is 
recovered and the perturbative solution is refined well. Then we show that this improvement 
method can apply for a ACDM model and improved the power spectrum of the density field. 

PACS numbers: 04.25.Nx, 95.30.Lz, 98.65.Dx 
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I. INTRODUCTION 



There are various structures in the universe which are gravitationally bounded, for ex- 
ample, galaxies, groups of galaxies, clusters of galaxies, voids, large-scale structure, and so 
on. These structures have evolved spontaneously from a primordial density fluctuation. 

The scenario for the growth of density perturbation is analyzed by several methods. When 
we do not consider the superhorizon scale or an extremely dense region like a supermassive 
black hole, the motion of the cosmological fluid can be described by Newtonian cosmology. 
Further, the Lagrangian description for the cosmological fluid can be usefully applied to the 
structure formation scenario. This description provides a relatively accurate model even in 
a quasi-linear regime. Zel'dovich jj| proposed a linear Lagrangian approximation for dust 
fluid. This approximation is called the Zel'dovich approximation (ZA) 



. After that, hig 
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rer-order approximation for the Lagrangian description was proposed 
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How accurate is the Lagrangian perturbation? To verify its accuracy, we often use simple 
models to compare exact and perturbative solutions. One of the simplest models is the 
"top-hat" spherical symmetric model, which has a constant density. For this model, we have 
obtained an exact solution. Therefore, to estimate validity in some approximated model, we 



often use the top-hat model. According to recent analyses of several symmetric models [18], 
the spherical symmetric model has a little difficulty accurately describing evolution with the 
Lagrangian perturbation. 



Munshi, Sahni, and Starobinsky [19( derived up to the third-order perturbative solution. 



20j obtained up to the fifth-order perturbative 



In addition to these, Sahni and Shandarin 

n 

solution. We [8j have derived up to the eleventh-order solution. 

If the density fluctuation is positive, the model will collapse. After a caustic formation 
at the center of the model, the equation of motion cannot describe the evolution. In past 
analyses with the Lagrangian perturbation, if we consider higher perturbative solutions, the 
approximation is improved all the more. On the other hand, if the density fluctuation is 
negative, the spherical void expands. In this case, ZA remains the best approximation to 
apply to the late-time evolution of voids. Especially if we stop expansion until an even order 
(2nd, 4th, 6th, • • ■) is reached, the perturbative solution describes the contraction of a void. 
From the viewpoint of the convergence of the series, we conjectured that the higher the order 



2 



of perturbation we consider, the worse the approximation becomes in late-time evolution. 
In this paper, we apply Shanks transformation, which accelerates the convergence of the 



sequence 



2l|. When we regard the perturbative solution as partial sums of infinite sequence, 



we must consider sequence convergence if description accuracy is to be discussed. Applying 
Shanks transformation, we can converge the sequence with a few terms. Therefore, we can 
improve the perturbative solutions. Here we consider the spherical void evolution using 
Shanks transformation. As a result, the higher-order perturbative solution recovers its ac- 
curacy and describes late-time evolution well. Because this method can apply not only to 
a spherical model but also in a more generic case, we can think that a new perturbative 
approach for Lagrangian description has been found. From comparison of Shanks transfor- 
mation and Pade approximation, which is another improvement method for the sequence, 
we show several merits of Shanks transformation for the perturbative approach. 

To generalize our approach, we apply Shanks transformation for a ACDM model. From 
the power spectrum of the density field, we show that Shanks transformation also recovers 
the accuracy of Lagrangian perturbation in a ACDM model. 

In Sec. mi we briefly show the evolution equation for a spherical symmetric model and 
derive exact solutions. Then we introduce Lagrangian perturbative solutions (Sec. IIIII) . 
In Sec. IIVt we describe Shanks transformation, the important method we have applied. 
Using this transformation, we indicate the accuracy of new perturbative solutions and show 
that Shanks transformation improves their solutions (Sec. |V|). For comparison, we also 
compute Pade approximation and show its behavior (Sec. IVT1) . In Sec. IVHl we consider the 
generic case and the evolution of a ACDM model with N-body simulation and Lagrangian 
approximations. Then we compare the power spectrum of the density field between the 
simulation and the approximations. Finally we offer our summary and conclusion (Sec. IVIIII) . 



II. EVOLUTION EQUATION AND EXACT SOLUTIONS 



We consider the "top-hat" spherical symmetric model, which has a constant density. In 
the E-dS Universe model, the equation of motion of a spherical shell is written as 



where x is a comoving radial coordinate and x$ = x(to) [19j]. Under the initial condition 
\5\ = a for a — > 0, Eq. ([T]) can be integrated. 



'dR\ _ n 3 
d^J ~ a \R~5 



(2) 



where R(9) = a(t)x/xo is a physical particle trajectory. The exact solution for the spherical 
collapse (Eq. (T5])) can be parameterized as follows: 



R+(0) = io(l-oosi 



a(9) 



sm i 



2/3 



(3) 
(4) 



Similarly, the exact solution for the expansion of a top-hat void (Eq. (T5])) can be parame- 
terized as follows: 



R46) = — (coshfl- 1) 



a(9) 



(sinh 0-0) 



-■ 2/3 



(5) 
(6) 



From these equations, we can obtain density fluctuation. 

*(x)=^V-l, 



x J 



(7) 



S + (x) 



9(6- sin Of 
2(1 - cos Of ' 
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-1, 
- 1 



(8) 
(9) 



where subscript + and 
spectively. 



2(cosh0- l) 3 

denote the case of spherical collapse and of void expansion, re- 



III. LAGRANGIAN PERTURBATION 



In the Lagrangian description, the inhomogeneity of mass distribution is described by 
the displacement from homogeneous distribution. The Lagrangian perturbative solution for 
spherical symmetric models in the E-dS Universe model is given by 



R±(t) = Ro 



l-Y^(±l) k C k a k 

k=l 



(10) 



TABLE I: The perturbative coefficients in Lagrangian description. 



k 
1 
2 
3 
4 
5 
6 
7 



10 
11 



1/3 
1/21 
23/1701 
1894/392931 
3293/1702701 
2418902/2896294401 
55964945/147711014451 
611605097/3430178002251 
4529700278678/52512595036460559 
2008868248800940/47103797747705121423 
29117328566723/1356899523596443827 



where C& are Lagrangian perturbative coefficients. The sign in coefficients corresponds to 
positive and negative density fluctuation, respectively. Substituting Eq. ( fill to ([2]), we 
derive the coefficients Ck- 

Mo-, Sa hni , and Star o binsky Q derived up t o th e third-order pe rtU r bati ve so.ution 
(Ci,C 2 ,C3). In addition to these, Sahni and Shandarin [20j obtained C 4 and C 5 . Further- 
more we (Si derived C*6, • • • , Cn. The coefficients Ck are shown in Table [fl 

The Lagrangian perturbation causes a serious problem. If the density fluctuation is 
positive, the spherical fluctuation collapses. The higher-order Lagrangian approximation 
gives accurate description (Fig. [p. On the other hand, if the density fluctuation is negative, 
the spherical void expands. In this case, ZA remains the best approximation to apply to the 
late-time evolution of voids. Especially, if we stop to improve until an even order (2nd, 4th, 
6th, • • ■), the perturbative solution describes the contraction of a void. In other words, ZA 
gives the best description for the late-time evolution of voids (Fig. [2]). 
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FIG. 1: The evolution of the spherical model for positive fluctuation. The thick solid line shows 
evolution by exact solution. The fine solid line, the dashed-dotted line, and the dashed line show 
evolution by first-, fifth-, and eleventh-order Lagrangian perturbation, respectively. In this case, 
the higher-order Lagrangian approximation gives an accurate description. In this figure, the scale 
factor is normalized by the time of a caustic formation in a first-order perturbation. 
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FIG. 2: The evolution of the spherical model for negative fluctuation. In this figure, the scale factor 
is normalized as in Fig. [TJ The thick solid line shows the evolution by the exact solution. The 
fine solid line, the dashed-dotted line, and the dashed line shows the evolution by first-, fifth-, and 
eleventh-order Lagrangian perturbation, respectively. In this case, the higher-order Lagrangian 
approximation deviates from the exact solution at a late time. In other words, ZA gives the best 
description for the late-time evolution of a void. 

IV. SHANKS TRANSFORMATION 

The Lagrangian perturbation causes a serious problem at late-time evolution for the 
spherical void model. How do we improve the description for spherical void evolution? 
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From the viewpoint of series, the higher order the perturbation we consider, the narrower the 
radius of the convergence becomes. How do we improve the convergence of the perturbation? 

For the contradiction in the Lagrangian approximation, we consider to improve the con- 
vergence rate of a sequence of partial sums. As a good way to speed up the convergence of a 
slowly converging series, Shanks transformation had been proposed [2l|. First we consider 
a simple example. Suppose the n-th term in the sequence takes this form: 

A n = A + aq n {\q\ < 1). (11) 

The sequence converges A n — > A as n —>■ oo. To obtain the limit of a sequence A, we solve 
algebraic equations with A n -i,A n , and A n+ \. 

A = ^ n + 1 ^ n ~ 1 ~ (12) 
A n +i + A n -\ — 2A n 

This formula is exact only if the sequence A n is described by the form in ffTTj) . For the 
generic case, we consider the nth term in the sequence takes the form: 

A n = A(n)+aq n , (13) 

where for large n, A(n) is a more slowly varying function of n than A n . Let us suppose that 
A(n) varies sufficiently slowly so that A(n — l),A(n), and A(n + 1) are all approximately 
equal. Then the above discussion motivates the nonlinear transformation 

S (A n) = A A ^ A i ;f; A . (M) 

This transformation is called Shanks transformation, creating a new sequence S(A n ) which 
often converges more rapidly than the old sequence A n . The sequence S 2 (A n ) = S[S(A n )] 
and S 3 (A n ) = S[S[S(A n )}} may be even more rapidly convergent. 

Damour, Jaranowski, and Schafer app lied Shanks transformation to the post-Newtonian 



approximation of general relativity 22j. They improved the analytical determination of 
various last stable orbits in circular general relativistic orbits of two point masses. In the 
next section, we apply the transformation to improve of the Lagrangian perturbation. 

V. IMPROVEMENT OF THE LAGRANGIAN PERTURB ATIVE SOLUTION 

Here we apply Shanks transformation for the Lagrangian perturbation. We consider 
the spherical void case and adopt the eleventh-order solution [8|. From the Lagrangian 
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FIG. 3: The evolution of the spherical model for negative fluctuation. In this figure, the scale 
factor is normalized as in Fig. [TJ We apply Shanks transformation to the Lagrangian perturbative 
solution and improve its accuracy. The thick solid line shows the evolution by the exact solution. 
The fine solid line shows the behavior of original eleventh-order perturbative solution. The dashed- 
dotted line, and the dashed line shows the evolution by once transformed and twice transformed 
perturbative solutions, respectively. After three times transformation, because the difference be- 
tween exact solution and perturbative solution becomes quite small. Therefore, if we also draw the 
curve of three times transformed perturbative solutions, we cannot find the difference on the figure. 
In this case, by application of Shanks transformation, the Lagrangian approximation improves its 
accuracy. In other words, we can describes late-time evolution of voids well. 

perturbation (Eq. (ITU]) ), we obtain a new solution via Shanks transformation. 



From the new sequence or perturbative solution R nj we can derive a more refined solution. 
Figure [3] shows the evolution of the spherical void using exact and the Lagrangian per- 
turbative solutions. We apply Shanks transformation once, twice, and three times for the 
perturbation. After transformation, the solution is refined and recovered its accuracy. In 
ordinary Lagrangian perturbation, we cannot improve the perturbative solution for late-time 
evolution. Using Shanks transformation, we can obtain a well-refined perturbative solution. 

We have shown one of the simplest cases. This improvement method is not limited to 
a special case. It can be applied in generic cases. Suppose the third-order perturbative 



Rn+lRn-l ~ R n 



(15) 




71 



(16) 



S 



solution takes this form: 



S = D + {t)S {1) + D+(t) 2 S {2) + D + (t) 3 S i3) 



(17) 



where D + is a linear growing factor. In the E-dS model, when we consider only the primordial 
growing mode in the longitudinal mode, the perturbative solution can be described by this 
form exactly 3,Q,Q- 



V 2 S (2) 



(v 2 s (1) ; 



c(l) 



V 2 ^ = -\det(S^) 



3 
10 

+ 21 



V 2 3 (1) V 2 3 (2) _ S (j) S (2) 



(18) 



(19) 



where V and subscript denote the Lagrangian spacial derivative. For other universe models, 
the perturbative solution can be approximated by this form in the matter-dominant era. 
Applying Shanks transformation (fill) , the perturbative solution is transformed to 

\2 



In Sec. IVIH we will treat a generic case with Shanks transformation. 



(20) 



VI. COMPARISON WITH PADE APPROXIMATION 

For a convergence of series, there are other methods. One of these methods is known 
as Pade approximation 21(. Pade approximation seems to be a generalization of Taylor 
expansion. For a given function /(£), Pade approximation is written as the ratio of two 
polynomials: 



f(t) 



(21) 



where otk and /?& are constant coefficients. Assume we already know the coefficient 7; (0 < 
I < M + N) of the Taylor expansion around x — 0. Then, 



M+N 

/(*)= E iit l + o(t 
1=0 



M+N+l\ 



(22) 



Comparing the coefficients (3%, and 7^, we determine and /5fc. 



a = 70 



(23) 



9 



N 

a k = E^mTfc-m (k = 1, ■ ■ ■ , N) , (24) 

m=l 

N 

PmlN-m+k = -lN+k (k = 1, • • • , N) . (25) 

m=l 

The advantage of Pade approximation is that even if we consider a same-order expansion, 
Pade approximation describes original function rather better than Taylor expansion does. 

Yoshisato, Matsubara, and Morikawa [3] have proposed an application of Pade ap- 
proximation 



or Eulerian perturbative solutions. Furthermore, Matsubara, Yoshisato, and 
Morikawa 24J have applied Pade approximation for the Lagrangian description. They also 
showed that Pade approximation can improve the Lagrangian perturbative solution. 

Here we apply Pade approximation to the spherical void. In Pade approximation, it is 
quite important to choose the numbers of terms M and N. When N is greatly different 
from M, the approximation is not improved well. Fig. 0] shows the evolution of the spherical 
void using exact and the Lagrangian perturbative solutions. We apply Pade approximation 
with several cases. Here we show the case of (M,N) = (1,10), (3, 8) and (5,6). It is 
important to choose the parameter M and N. When N is greatly different from M, not 
only the approximation is not improved well, but also the solution will diverges. For late-time 
evolution, although the original eleven-order perturbative solution converges to 5 — ► —1, the 
Pade approximation with (M,N) = (1,10), (3, 8) diverges. On the other hand, when we 
choose (M, N) = (5, 6), the perturbative solution can approximate the exact solution at late 
time. We can improve the perturbative solution with Pade approximation, too. 

VII. GENERIC CASE: THE ACDM MODEL 

We showed the evolution of a homogeneous spherical void and noted the improvement for 
Lagrangian approximation. Because we know the exact solutions for the spherical collapse 
and void evolution, we do not know whether or not the Shanks transformation is useful for 
the generic case. Therefore we must apply the transformation for generic models. 

Here we consider a ACDM (Low density Cold Dark Matter) model. The cosmological 
parameter at the present time (z = 0, here we define a = 1 at the present time) is given by 



a WMAP 3rd-year result 



25|: 



tt m = 0.28, (26) 
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FIG. 4: The evolution of the spherical model for negative fluctuation. In this figure, the scale factor 
is normalized as in Fig. [TJ We apply Pade approximation to the Lagrangian perturbative solution 
and improve its accuracy. The thick solid line shows the evolution by the exact solution. The fine 
solid line shows the behavior of original eleventh-order perturbative solution. The dashed-dotted 
line, and the dashed line shows the evolution by the case of (M, N) = (1, 10) and (3, 8), respectively. 
When N is greatly different from M, the approximation is not improved well. When we choose 
(M,N) = (5,6), the difference between exact solution and perturbative solution becomes quite 
small. Using Pade approximation, we can improve the perturbative solution well. 



&DE 



0.72, 

74 [km/s/Mpc] 
0.74. 



(27) 
(28) 
(29) 



The Gaussian density field is generated by COSMICS [26(. We set up the initial condition 
at decoupling time (a = 10~ 3 ). The initial peculiar velocity and the density fluctuation are 
adjusted by the growing solution of ZA. 

For time evolution, we consider Lagrangian third-order approximation, Shanks transfor- 
mation, and N-body simulation. For computation of the Lagrangian perturbations, we set 
the parameters as follows: 

Number of grids : N = 128 3 , 

Box size : L = 128/i _1 Mpc (at a = 1) . 



The N-body simulation is applied by a particle-particle particle- mesh (P 3 M) method 27 ] 
whose code was originally written by Bertschinger. 
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For N-body simulations, we set the parameters as follows: 

Number of particles : iV = 128 3 , 

Box size : L = 128/i _1 Mpc (at a = 1) , 
Softening length : e = 50/i _1 kpc (at a = 1) . 

Then we impose a periodic boundary condition. 

The Lagrangian approximation in ACDM is expanded as 

S = h!{t)SW + h 2 (t)S^ + h 3 (t)S® , (30) 

where h n (t) is the growing factor for nth-order approximation. The spacial parts are given by 
Eqs.f llSp and f)19p . The growing factors are derived with a numerical method Applying 
Shanks transformation (Eq. (114jl ). the perturbative solution is transformed to 

§ = hlSil) + h 3 S<?> - hsw ■ (31) 

In order to avoid the divergence of the density fluctuation, we need to consider a smoothed 
density field over the scale R. This density field is related to the unsmoothed density field 
p(x) as 

p(x;R) = J d 3 yW(\x-y\;R)p(y) 
d 3 k T ~ T 

W) 

where W denotes the window function and W and p represent the Fourier transforms of the 
corresponding quantities. In this paper, we adopt the top hat window function, 

^ = 3(sinx-xcosx) 

X 3 

Then, the density fluctuation S(x; R) at the position x smoothed over the scale R can be 
constructed in terms of p(x; R). 

Here we choose the smoothing scale R = l/i^Mpc. Then we calculate the power spectrum 
of density fields. In order to obtain the power spectrum, we generate 50 samples for the 
primordial density fluctuations. Then we pick up snapshots at z = 5, 4, 3, 2, 1, 0. 

Figure [5] shows the power spectrum of the density field in N-body simulation. During 
evolution, because a strongly nonlinear region promotes the growth of the fluctuation, the 
large-fc components in the power spectrum grows remarkably. 
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/ ^W(kR)p(k)e- ih - x , (32) 
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FIG. 5: The power spectrum of the density field in N-body simulation. Because of a nonlinear 
effect in small structures, the large-/c components in the power spectrum grow remarkably. 
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FIG. 6: The power spectrum of the density field in Lagrangian third-order approximation. Because 
of the shell-crossing, the spectrum distorted from the large-A; components to the small-A; components 
gradually. At z = 0, the spectrum is really differs greatly that in N-body simulation. 

On the other hand, the power spectrum in the Lagrangian third-order spectrum sinks 
(Fig. [6]). At the high-z era, shell-crossing occurs in small structures. After shell-crossing, 
the cluster it forms spreads eternally. Then the negative influence of shell- crossing affects 
the large-scale structure. Therefore the spectrum sinks from the small scale gradually. At 
z = 0, the spectrum differs greatly from that in N-body simulation. 

When we apply Shanks transformation, although a strongly nonlinear effect in the large-/c 
components is not realized, the distortion of the spectrum is well improved (Fig. [7]). Even 
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FIG. 7: The power spectrum of the density field in Shanks transformation. Because of the im- 
provement of the nonlinear effect, the distortion of the spectrum is improved well. 

if we consider a low-z era, Shanks transformation can realize the spectrum, except for the 
small scale, well. 

VIII. SUMMARY 

We have discussed the evolution of the spherical void in the framework of the Lagrangian 
perturbation. Using ordinary Lagrangian perturbation, the higher-order Lagrangian approx- 
imation deviates from the exact solution at late time. In other words, ZA gives the best 
description for the late-time evolution of voids. Then we generalized Shanks transformation 
for the Lagrangian perturbation, i.e., we analyzed time evolution for a ACDM model. 

We apply Shanks transformation, which accelerates the convergence of series for the 
Lagrangian perturbation. Although the transformation in valid within a convergent radius 
of the series, the transformation creates a new sequence which often converges more rapidly 
than the old sequence. In the spherical void model, the transformation is valid for a long 
time. Then we can improve the accuracy of the Lagrangian description. Using this method, 
we can solve the problem whereby the higher order the perturbation we consider, the worse 
the approximation becomes in late-time evolution. 

In this paper, we also compare the accuracy of our improved methods, between Shanks 
transformation and Pade approximation, that is. In the comparison, we found that Shanks 
transformation has several merits: both Shanks transformation and Pade approximation 
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can be derived with algebraic procedures. In Pade approximation, although we can solve 
algebraic equations and find unique solutions, the equation is quite complicated. Further- 
more, according to our analyses, if the difference between two parameters M and N in Pade 
approximation is large, the perturbative solution will diverge. On the other hand, Shanks 
transformation does not diverge. 

However, in several points, Shanks transformation shows its weakness. To apply Pade 
approximation, we need to obtain second-order perturbation. On the other hand, when we 
consider the Shanks transformation, we must obtain at least third-order perturbation. Using 
the Shanks transformation, we obtain a new perturbative solution R n from R n -i, R ni R n +i- 
Then, to repeatedly apply Shanks transformation, we require R n -2, • ■ • , R n +2- m general, 
when we apply n times transformation, at least we must know 2n + 1-th order perturbative 
solutions. To improve the perturbation well, we must repeat the transformation several 
times. 

From the viewpoint of algebraic procedures, Shanks transformation has an advantage. In 
Pade approximation, we must solve nonlinear simultaneous equations. Then the solution is 
extremely complicated in a higher-order case. For example, when we improve an eleventh- 
order solution with (M,N) = (5,6), we derived about fifty-digit coefficients. 

For a ACDM model, Shanks transformation recovers accuracy for the description of the 
density field, too. We showed that the power spectrum in Shanks transformation becomes 
better than that of ordinary Lagrangian approximations. Even if we consider the spectrum 
at z = 0, the Shanks transformation can describe quite well, except for the small scale. 

However, in a generic case, the critical problem in Shanks transformation appeared. 
When we continue to apply Shanks transformation, the divergence of the perturbation oc- 
curs. For example, when we apply the transformation with third-order approximation, the 
perturbation is written as Eq. (1511) . If the higher-order perturbation is smaller than the 
lower-order perturbation, we consider that the perturbative method is valid. From Eq. (l3~Tl) . 
when second-order and third-order perturbation become 

h 2 S {2) ~ h 3 S {3) , (34) 

the perturbation diverges. In other words, because the Lagrangian perturbation is excluded 
from the convergence series during evolution, the perturbative method becomes invalid. 
Therefore when we apply Shanks transformation, we notice the validity of the perturbative 
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expansion. In the spherical collapse case, the density fluctuation diverges before divergence 
of the Lagrangian perturbation. After the shell-crossing, the divergence of the perturbation 
is settled. Then the perturbation describes the displacement of the fluid. 

Even if we improve the Lagrangian perturbation, we cannot avoid the problem of shell- 
crossing. Because of shell-crossing, the approximation becomes worse at late time. To solve 
this problem, Scoccimarro and Sheth proposed an extrapolation method j^]. First, they 
considered structure formation with second-order Lagrangian approximation. After that, 
they extrapolated the density distribution in the high-density region from NFW profile. 
This hybrid method realizes the density distribution well. Although it is very complicated, 
this method is effective for higher-order perturbation. We can then expect that this hybrid 
method will realize the density distribution with high precision by combining it with third- 
order Lagrangian approximation or Shanks transformation. 
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